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1.  INTRODUCTION  AMD  SUMMARY 


1.1  INTRODUCTION 


Radomes  can  limit  the  performance  of  airborne  radar  and  guidance  systems  by  reduc¬ 
ing  guidance  accuracy  and  radar  range.  This  limitation  persists  despite  continuing 
radome  development  because  newer  missiles  fly  foster  and  higher,  and  because  newer 
threats  require  greater  frequency  bandwidths.  In  addition,  the  need  for  maximum  ra¬ 
dar  range  drives  the  antenna  diameter  to  foe  maximum  that  fits  within  the  radome; 
consequently,  the  antenna  receives  waves  that  are  guided  by  the  radome.  These  guid¬ 
ed  waves  produce  bo  resight  errors.  Our  experiments  have  shown  that  the  guided 
wave  induced  errors  require  radome  correction  techniques  that  differ  from  the  usual 
method  of  small,  local  thickness  variations. 

The  electromagnetic  design  of  radomes  and  foe  correction  of  radome  boresight  error 
usually  requires  a  long  development  period  that  includes  much  empirical  work.  Al¬ 
though  analytical  methods  have  been,  and  k.*e  being,  developed,  most  methods  are  ap¬ 
proximate.  Unfortunately  systematic  error  analyses  are  limited  because  radomes  in¬ 
volve  many  factors  and  because  foe  enclosed  antenna  influences  performance. 

Because  of  these  limitations  on  radome  design  methods,  interest  exists  in  electronic 
compensation  methods.  Radome  boresight  error  would  be  measured  and  stored  in  an 
on-board,  digital  memory,  and  these  data  would  be  applied  as  corrections  during  mis¬ 
sile  flight.  This  approach  is  attractive,  but  it  has  limitations.  Obviously,  adequate 
storage  is  required.  The  memory  must  be  large  since  boresight  error  depends  on 
glmbal  angle,  frequency,  and  wave  polarization.  Existing  memories  and  processors 
seem  capable  of  handling  these  variables  within  size  limitations  of  missiles,  but  fre¬ 
quency  and  polarization  sensors  would  be  required.  Cost  impacts  need  study. 

A  potentially  more  serious  problem  is  foe  temperature  dependence  of  radome  bore¬ 
sight  error.  This  dependence  involves  many  possible  trajectories  and  the  missile  po¬ 
sition  within  each  trajectory.  However,  no  accurate  facility  exists  for  measuring 
boresight  error  of  heated  radomes;  this  is  an  obstacle  to  electronic  compensation. 

An  alternative  to  electronic  compensation  is  designing  the  radome  for  elevated  tem¬ 
peratures.  The  radome  would  be  somewhat  thinner  than  for  best  room  temperature 
performance.  However,  this  alternative  still  requires  high  temperature  measure¬ 
ments  of  boresight  error  for  its  evaluation,  and  it  depends  on  accurate  design  methods. 

Accurate  radome  design  methods  require  accurate  descriptions  of  wave  propagation 
through  foe  radome.  This  description  must  Include  the  effects  of  frequency,  polariza¬ 
tion  and  foe  direction  of  the  incident  wave.  The  radii  of  curvature  of  the  radome  are 


{ 


also  a  significant  factor  because  they  influence  the  wave  types  excited  fay  an  incident 
wave  and  the  accuracy  of  approximate  descriptions.  Typical  radii  range  from  less 
than  one  to  several  wavelengths. 

This  report  describes  a  continuing  study  of  wave  propagation  through  dielectric  shells. 
Emphasis  this  year  was  on  axially  symmetric  shells  such  as  hollow  cylinders.  Last 
year  we  considered  the  hollow  wedge,  which  is  a  two-dimensional  radome,  and  die  flat 
slab,  which  is  a  building  block  of  a  wedge  and  a  local  approximation  to  a  shell.  ** 
Quantitative  data  were  given  on  phase  shifts  caused  by  refraction  and  plane  wave  prop¬ 
agation  through  the  dielectrics;  in  addition  some  guided  wave  properties  were  deter¬ 
mined.  The  analytical  method  used  was  the  moment  method2,  . 

1.2  SUMMARY 

**  I  i 

Sections  describes  a  technique  estimating  the  amplitude  and  phase  of  guided  waves 
excited  on  finite  dielectric  slabs  by  incident  plane  waves.  The  procedure  is  based  on 
the  moment  method,  but  it  extends  that  method  by  postulating  an  explicit  form  for  the 
guicted  waves  and  solving  for  the  guided  wave  amplitudes.  The  method  solves  for  the 
two  guided  waves  that  travel  in  opposite  directions  parallel  to  the  slab.  Experiment 
verified  the  technique. 

Section#  describes  nearfield  scattering  by  a  finite,  solid,  dielectric  cylinder.  Inci¬ 
dence  was  broadside.  The  procedure  extends  the  theory  originally  developed  by 
Richmond.  It  identifies  the  effects  from  the  ends  of  the  cylinder.  Experiment  verified 
the  procedure  for  polarization  parallel  or  perpendicular  to  the  cylinder. 

. 

Section  4  describes  a  new  technique  for  nearfield  scattering  by  a  hollow  dielectric 
cylinder.  The  moment  method  was  used  but  the  cells  were  finite  angular  sectors. 

This  is  a  new  shape.  Singularities  of  integrands  were  integrable.  Experiment  verified 
our  approximate  theory,  especially  when  a  correction  for  the  pattern  of  the  probe  an¬ 
tenna  was  utilized. 

■V- 


♦References  are  given  in  Section  5. 
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FINITE  DIELECTRIC  SLAB:  GUIDED  WAVE  AMPLITUDE 


2.1  INTRODUCTION 

This  section  describes  an  application  of  the  moment  method  to  determining  the  guided 
wave  magnitude  on  a  finite  dielectric  slab.  Reference  1  presented  a  moment  method 
analysis  of  flat  slabs,  and  it  gave  experimental  checks.  The  method  will  be  used  in 
this  report,  but  only  a  brief  summary  of  the  theory  is  included. 

2.2  THEORY 

Let  us  represent  the  field  within  the  slab  as  a  sum  of  four  waves,  which,  in  the  co¬ 
ordinate  system  of  Figure  2-1,  are  as  follows: 

1)  A  plane  wave  transmitted  into  the  slab:  TeikxX, 

2)  A  reflected  plane  wave:  Re^x*, 

3)  A  slab-guided  wave:  Ae^X, 

4)  A  reflected  slab-guided  wave:  Be  i^tSx, 

At  the  slab  midplane ,  Y  equals  zero;  thus  the  field  is 

ik  x  lk  x  -ik  x 

E(y  =  0)  =  Se  X  +Ae  g  +Be  8  ,  (2-1) 
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Figure  2-1.  Co-ordinates  for  a  Slab 
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where  by  assumption  S  is  T  +R,  and  T  and  R  are  respectively  the  complex-valued 
transmittance  and  reflectance  of  a  flat  dielectric  sheet  for  plane  waves.  In  addition, 
kx  is  (2ir/X)  sin  a ,  with  X  the  wavelength  and  a  the  Incidence  angle,  and  kg  the  propa¬ 
gation  constant  of  a  wave  guided  by  an  infinite  slab.^  We  seek  A  and  B. 

The  total  field  in  a  slab  can  be  accurately  determined  with  the  moment  method.  Be¬ 
cause  that  theory,  extensive  calculations,  and  measurements  were  given  in  Reference 
1  through  3,  we  only  summarize  the  method  here.  The  moment  method  starts  by  di¬ 
viding  a  slab  into  a  set  of  N  parallel  cylinders  and  assuming  the  field  constant  within 
any  cell.  An  integral  equation  for  scattering  relates  the  unknown  total  field  E^,  the 
incident  field  E1,  and  the  scattered  field  Es,  so  that 

ET  -  (ik2/4)  ( k  - 1)  SS  E-  (x»,  y»)  H  (1)  (kp)  dx*  dy*  =  E1 ,  (2-2) 

o  i 

where  H0^  is  a  Hankel  function,  k  is  dielectric  constant,  and  p  is  distance  between 
observation  and  integration  points.  Equation  2-2  is  evaluated  at  the  center  of  each 
cell,  with  the  integral  extending  over  all  cells,  to  produce  N  simultaneous,  linear  al¬ 
gebraic  equations.  The  result  is  that 


I 


n=l 


mn  n 


(2-3) 


for  1  2S  m  s  N. 

To  determine  A  and  B,  combine  Equations  2-1  and  2-3  and  use  Snell's  law.  Thus, 


SLC  e 
mn 
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x  n 
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mn 


ik  x 
g  n 


+  B£c  e 
mn 


■ik  x 

g  n 


=  e 


ik  x 
x  m 


(2-4) 


for  m  equal  1  through  N.  Now  N  can  be  large.  Its  value  depends  on  foe  size  of  foe 
slab  because  cylinder  radius  is  limited  by  foe  assumption  that  foe  field  is  constant. 
Although  N  is  large,  we  have  only  two  unknowns  (A  and  B)  because  these  guided  wave 
amplitudes  were  assumed  constant.  Actually  A  and  B  could  depend  on  x.  In  this  case 
we  would  have  N  values  of  A  and  N  values  of  B.  However,  if  A  and  B  vary  slowly  with 
x,  it  seems  reasonable  to  consider  successive  pairs  of  equations  in  Equation  2-4  and 
to  solve  each  pair  for  a  value  of  A  and  a  value  of  B.  The  solution  is  of  course  approx¬ 
imate  because  higher  order  modes  can  exist,  and  even  foe  dominant  slab  mode  may 
not  be  fully  excited  near  an  edge. 


To  evaluate  A  and  B  requires  first  evaluating  foe  co-efficients  Cmn  and  then  solving 
successive  pairs  from  Equation  2-4.  For  m  and  m +1  we  obtain 
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A  =  [r  -3*  r  (s’  ,)_1]  Is  -s’  s  (s’  _)’1) , 
m  m  m+1  m+1  m  m  m+1  m+1 


where 


s  =  S  c 
m  L  mn 

e  g*. 

(2-6) 

s  =Sc 

m+1  m+1 

ik  x 

e  gn 
,n® 

(2-7) 

s'  =  >C 

m  '  mn 

-ik  x 
e  g“, 

(2-8) 

-ik  x 

s'  =Yc  e  gn, 

m+1  Z  m+l,n 

(2-9) 
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x  m  _  \ 
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L  mn 

(2-10) 

ik  x  , 
x  m+1 

r  =  e  -  S  ) 

m+1  Z 

_  Ik  x 

_  x  n 

►C  ,  e 
-  m+1,  n 

(2-11) 

In  addition. 


B  =  <rm+l  '  A  W  <sm+l> 


(2-12) 


2.3  COMPUTATION 


To  test  the  approximate  solution  we  analyzed  a  slab,  with  dielectric  constant  2.6, 
thickness  0.25”,  and  length  4”.  Experimental  verification  for  the  slab  was  given  in 
Reference  3.  Although  the  existence  of  guided  waves  was  established  and  magnitudes 
were  estimated,  no  explicit  values  of  A  and  B  were  given.  Frequency  was  9.375  GHz 
and  polarization  was  perpendicular.  Figure  2-2  shows  computed  values  of  A  and  B. 
We  see  that  computed  values  of  A  increase  with  x  for  grazing  incidence  *  0  (or  a 
=  90*),  but  A  decreases  with  x  for  $  -  22. 5\  For  the  other  values  of  $  ,  the  values 
oscillate.  1 
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from  Equations  2-5  and  2-12) 


2.4  MEASUREMENT 


The  field  near  the  slab  was  measured  with  a  microwave  interferometer  that  utilized  a 
dipole  probe.  Figure  2-3  shows  an  example  of  measured  data  for  4>l  -  30*;  the  probe 
was  spaced  0. 1"  from  the  slab.  We  estimate  guided  wave  amplitude  A  by  first 
smoothing  the  curves  to  delete  the  reflected  guided  wave,  amplitude  B.  We  obtain  two 
values  of  A,  at  minima  and  maxima  on  the  condition  that 


E  =  T  +Ae~vd  (2-13) 

P 

where  Tp  is  the  plane  wave  power  transmittance,  v  is  the  decay  co-efficient  for  the 
field  outside  the  slab  and  d  is  perpendicular  distance.  We  utilize  the  extrema  because 
at  Ihese  features  the  two  terms  on  the  right  hand  side  of  Equation  2-13  are  either  in  or 
out  of  phase.  We  then  subtract  Tp  (which  is  easily  calculated)  from  |  E  [ ,  which  is 
known  from  the  measured  graph  of  power  J  E  | 2 .  After  subtraction,  the  difference  E 
-  Tp  is  corrected  by  multiplying  it  by  ev<*.  We  calculated  v  as  1. 15  ir/inch  from  the 
theory  of  Reference  1.  Figure  2-4  shows  the  values  of  A  determined  from  the  meas¬ 
ured  data  for  five  incidence  angles. 

2.5  DISCUSSION 

In  Figure  2-4  the  experimentally  determined  values  of  A  fall  within  the  envelope  of  the 
values  of  A  computed  by  solving  pairs  of  equations  from  Equations  2-4.  This  agree¬ 
ment  verifies  the  approach  to  solving  the  guided  wave  excitation  problem  for  a  slab  of 
finite  length.  The  solution  shows  that  A  varies  with  position  even  for  a  single  inci¬ 
dence  angle.  Our  verification  with  experiment  shows  the  model  and  technique  are 
reasonable.  We  point  out  that  we  did  not  compare  computed  values  of  A  at  specific  lo¬ 
cations  (values  of  x)  with  the  experimental  values  at  the  same  locations.  However, 
the  good  agreement  shown  in  Figure  2-4  substantiates  the  model  in  Equation  2-1  *  and 
the  approximate  pairwise  solution  of  Equations  2-4. 

In  summary,  we  have  utilized  the  accurate  moment  method  as  a  step  in  a  new  proce¬ 
dure  that  estimates  explicitly  the  magnitudes  of  guided  waves  on  dielectric  slabs.  Al¬ 
though  the  moment  method  gave  the  total  field,  it  did  not  identify  constituent  waves  of 
distinct  types.  This  identification  may  become  useful  in  correcting  ra domes  to  reduce 
boreslght  error  because  the  common  radome  practice  of  changing  wall  thickness  to 
correct  refractive  errors  seems  ineffective  for  errors  caused  by  guided  waves.  It 
seems  possible  that  the  method  can  be  extended  to  replace  A  and  B  with  functions  of  x 
while  retaining  the  functional  behavior  given  by  Equation  2-1.  Another  possibility  is 
to  utilize  N/2  values  of  An  and  N/2  values  of  Bq,  on  the  assumption  that  An  is  identi¬ 
cal  for  an  adjacent  pair  of  cells  and  similarly  for  Bq.  We  did  not  pursue  a  more 


*This  model  was  described  in  Reference  3. 


Figure  2-3.  Measured  Intensify  Near  Slab  (Polarization  Perpendicular,  Dipole  Probe 
0,1"  from  Slab,  Incidence  Angle  60°;  =  30°).  The  dashed  extension  of  that 

measured  with  slab  is  a  smoothed  version.  Arrows  show  the  transmittance 
at  peak  intensity  and  (he  minimum  value  of  intensity. 
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INCIDENCE  ANGLE  (OEG) 


Figure  2-4.  Guided  Wave  Amplitude  (Computed  from  equation  2-5  (> ); 
n  denotes  the  edge  nearest  the  source,  f  denotes  edge  further 
from  source,  measured  at  maxima  (+),  at  minimum  (-); 
data  corrected  by  exp  v  d 


detailed  evaluation  of  the  accuracy  of  this  approach  because  the  results  in  Figure  2-4 
are  reasonable  and  because  the  effects  of  curvature  seemed  more  significant.  The 
next  sections  introduce  curvature;  of  course,  they  also  consider  finite  objects. 
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NEAR  FIELD  SCATTERING  BY  FINITE,  SOLID  CYLINDER 


3.1  INTRODUCTION 

This  section  describes  nearfield  scattering  by  a  finite  length,  solid,  dielectric  cylin¬ 
der.  The  solid  cylinder  was  studied  because  it  may  be  useful  as  an  element  for  de¬ 
composing  a  hollow  cylinder  or  an  axially  symmetric  shell.  The  analysis  extends  that 
of  Richmond  for  infinitely  long  cylinders. 3 

3.2  THEORY 

The  field  scattered  by  a  dielectric  object  can  be  expressed  as  a  volume  integral  over 
polarization  currents.  The  scattered  field  is 

I®  <rf)  =  a.x  +  V  Lj)/4v  (3-1) 

where 

2  f  T 

Lj  =  *  I  (K-l)gl 

Lj  =  J gv  .  t(K-l)ET]  dv, 

and  Et  is  the  total  field,  g  is  r-1  expikr,  and  r  is  the  distance  between  integration  and 
observation  points.  The  prime  superscript  (on  the  gradient  operator  and  r’)  denotes 
the  observation  point. 

T 

The  moment  method  determines  E  ,  with  the  integral  in  Equation  3-1  and  the  relation 
that 


(3-2) 

(3-3) 


T  s  I 

E  -E  =  r  (3-4) 

where  E*  is  the  known,  Incident  field.  Equation  3-4  is  an  integral  equation.  It  is  re¬ 
duced  to  a  set  of  simultaneous  algebraic  equations  by  subdividing  the  scatterer  into 
small  cells  for  the  purpose  of  evaluating  Equations  3-3  and  3-4.  The  subdivision  jus¬ 
tifies  assuming  the  field  constant  within  each  cell;  fills  assumption  simplifies  the  inte¬ 
gration.  This  subdivision  leads  to  a  set  of  simultaneous  equations  from  Equation  3-1. 


The  following  sections  describe  the  theory  for  calculating  E  .  To  simplify  notation  let 

E8  (rf)  *  E8  +  E8  (3-5) 


where 


sA  -  V4'  • 


E  =  VI  Jit 

S  a 


3.  2. 1  EVALUATION  OF  E  J 

The  integral  Ij  is  expressed  as  a  sum  of  integrals  over  discs  like  those  in  Figure  3-1. 
Assuming  k  constant  throughout  the  cylinder  and  that  the  field  is  constant  within  any 
one  cell,  we  have  from  Equation  3-1  that 


Figure  3-1.  Subdivision  of  Cylinder  Into  Discs 


where  En  is  the  total  electric  field  in  the  disc  labelled  n.  Equation  3*1  leads  to  N  si 
multaneous  algebraic  equations  for  m  *  1,  .  .  .  N. 

To  evaluate  the  integrals  in  Equation  3-1  consider 
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2  2  e 
p  +  p’  -2  pp*  cos  (<p  -  $')  +  ( z  -  z’)‘ 


,1/2 


(3-9) 


Because  the  field  is  assumed  constant  within  a  disc,  the  field  is  evaluated  at  p'm  equal 
zero;  the  <f>  dependence  is  removed  so  the  integral  over  $  given  2rr.  We  have  then 
that 


V'  T 

+dz 
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r  2  r  21X^2 
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n  J 

r  2 

2-f1/2 

LP  4  (Zn-Z^)  J  dv 

where  a  is  the  cylinder  radius  and  each  disc  has  thickness  26z.  The  integral  over  p 
can  be  done  with  the  substitution  that 


u 


to  obtain 


*  -i2ff  k  (k  - 1)^>  EQ 

~z  +dz 
n 

r  2  -1V2 

exp  ik  [a  +(2-^)]  -  exp  ik  (z-z^) 

w 

z  -dz 

l 
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There  are  two  cases.  First,  z  *  z’  ;  in  this  case  expand  the  exponent 

n  m 


(3-11) 


in  a  binomial  series.  The  resulting  functions  have  the  form  x"1  sinx  where  x  is  a  lin¬ 
ear  in  fiz.  For  small  5z  x_1  sin  x  is  approximately  one  so  that 

N  1/2 

Lj  =  -i4ir  (k-1)  kfiz  ^  jexpik  ja2  +  -  expik  J  . 

n=l  (3-12) 

In  the  second  case  zn  =  z^;  z  -  z{q  is  replaced  by  die  variable  £  which  ranges  be¬ 
tween  5z  and  -  6z.  Again  expanding  the  exponential  in  Equation  3-11,  we  obtain  Fres¬ 
nel  integrals  which  we  approximate  by  a  small  argument  expansion,  to  obtain 


I,  =  4 ir  (jc  -  1)  (ka)  (kdz)  E1  . 

~m 


(3-13) 


Note  that  I_i  has  the  same  form  for  E  polarized  in  the  y  direction  and  in  the  z 
direction. 

The  scattered  field  is  found  from  Equations  3-12  and  3-13,  and  the  definition  in  Equa¬ 
tion  3-6. 

3. 2. 2  Eg  FOR  THE  FIELD  PARALLEL  TO  THE  CYLINDER  AXIS 
In  this  case,  E  has  only  a  z  component  so  Equation  3-3  gives 

'2  '/*f,"‘-1>EzdV-  (3-M> 

Let  the  length  of  the  cylinder  be  2L.  Then  (k  - 1)  is  discontinuous  with  jumps  at  z 
=  ±L.  We  write  it  as  (k-1)  [H  (z  +  L)  -  H  (z  -  L)] ,  where  H  is  the  Heaviside  unit  func¬ 
tion.  The  partial  differentiation  leads  to  Dirac  delta  functions.  We  assume  Ez  varies 
slowly.  Again  decompose  the  integral  into  a  sum  of  integrals  over  discs  like  in  Fig¬ 
ure  3-1,  with  the  field  constant  in  each  disc,  and  integrate  over  z,  to  obtain 


where 


(K-1) 


-1 


*2  * 


E  I  (-L)  -  E  I  (L) 
zl  V  '  zN  2l  ' 


(3-15) 


(exp  ikr_)  pdpd<l>  , 


(3-16) 
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and 


na  -1 

r+  (expikr+)  pdpdtp 
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2  2  2 
r  =  [p  +  p’  -  2pp'  cos  (4>  -  $')  +  (L  ±  z)  ] 


1/2 


(3-17) 


(3-18) 


To  evaluate  the  integrals  in  Equation  3-9,  we  expand  the  arguments  of  the  exponentials 
to  second  order.  Note  that  we  cannot  set  p'  to  zero  (as  in  the  previous  section)  be¬ 
cause  of  the  gradient  operator.  After  doing  the  integrals  and  the  gradient  operation 
we  obtain 


1/2  (fc-1) 


ikL+ 

e 


(3-19) 


and  die  scattered  field  from  I  (-  L)  is 


1/2  (K-l) 


I  2 

a  + 


ikL 

e 


(3-20) 


where  L+  is  (L  +  z’)  and  L  is  (L  -  z’). 

The  result  of  this  subsection  is  that  the  integrations  lead  to  two  terms  for  the  scattered 
field  from  Equation  3-13  and  two  terms  from  Equation  3-14.  These  terms  arise  from 
the  6  functions  associated  with  the  ends  of  the  cylinder.  Recall  dEz/dz  was  omitted. 

3.  2.  3  e|  FOR  THE  FIELD  PERPENDICULAR  TO  THE  CYLINDER  AXIS 


In  diis  case  the  field  has  p  and  <p  components  so 


?  ■  <3-21> 

As  before ,  we  assume  that  file  field  depends  only  on  z;  in  addition  (ft  - 1)  is  independ¬ 
ent  of  <p .  Therefore 
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Is  -  /g  £<«-!>£„<*.  (3-22) 

op  p 

Let  tiie  function  (k  - 1)  be  (k  - 1)  [1  -  H  (a)]  where  now  k  is  the  value  in  the  dielectric 
and  H  (a)  is  the  Heaviside  function.  Assume  the  field  has  only  a  y  component  so  Ep  is 
Ey  sin4>.  By  decomposing  the  integration  into  the  discs  of  Figure  3-1  we  have 

I2  =  -(K-llE^/gsin#  6(a)  dv.  (3-23) 

To  do  the  integral  first  integrate  on  p  and  then  introduce  the  variable  <f>  -  <p'  =  9.  We 
obtain 


\  -  -<«-l)as 


z  +6z  2ir-4>* 

r  i 

Z  -6zJ -4b* 


(cos <#>'  sin0  +  sin <p*  cos  0)  g(a)  dOdz 


(3-24) 


where  g(a)  is  g  in  Equation  3-3  with 


22  2  1/2 

r  *  r(a)  =  [a  +  p*  -  2ap*  cos  9  +  (z  -  z’)  ]  .  (3-25) 

The  first  term  on  the  right  side  gives  an  Integral  equal  zero  because  the  integral  is  an 
odd  function  of  6.  To  evaluate  the  second  term  expand  r  (a)  in  a  binomial  series 

r(a)  »  d  -  ad  *  p*  cos0  +  £d  1  (3-26) 

r  2  o-JL/2 

where  d  is  [(zn  -  z^)*  +  a  j  ,  and  is  (zn  -  z^).  Recall  zn  is  at  the  center  of  the 
disc  for  which  the  integral  is  being  done  and  z^  is  the  observation  point,  at  the  center 
of  one  of  the  discs.  The  Integration  yields 

I„  *  2i  it  (k-1)  sine  (k5z)  p*  sin^f  2  E  d  1  expikd  (3-27) 

2  yn 

where  sincx  indicates  x  1  sinx.  The  scattered  field  follows  from  Ihe  gradient  opera¬ 
tion  in  Equation  3-1;  the  y  component  is 

(ES)y  =  ~(x-l)  ka2  3z  sine  (kflz)  2  E^  d  1  expikd  (3-28) 
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We  see  that  the  scattered  field  from  the  gradient  of  the  scalar  potential  contains  con* 
tributions  from  all  discs.  This  result  contrasts  with  the  case  for  the  field  parallel  to 
the  cylinder.  In  that  case  there  were  only  contributions  from  the  ends  for  the  gradient 
of  the  scalar  potential. 

3.2.4  ALGEBRAIC  EQUATIONS  FOR  INTERNAL  FIELDS 
The  algebraic  equations  result  from  the  definition 


T  S  I 
E  -  E  =  E 


For  the  electric  field  parallel  to  the  cylinder  we  have  for  m  =  1, 

N 


V  C  E  =  E 
mn  zn  m 


where 


n=l 


C11  =  °NN  =  1  "  P  ‘  B  ’ 


cin  “  iyp  hn~A ' 
cni  =  1  n  Su  "  A  5 

for  Ihe  other  diagonal  matrix  elements, 


and  for  Ihe  remaining  elements 


C  =  1  -  p  , 
nn 


C  =  ly  13  f 
mn  mn 


where  y  =  k-1,  0  =  k6z,  a  =  ka,  p  =  a/9y,  and 


(3-29) 
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+  (zn 


-  expik  (z^ 
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mn 


In  addition 


:u 
'  ic 


L.n 


/r-.i 


A  »  (y  a2/4)  {ik  -  (2L  -  fiz)”1]  (2L  -  Szf1  expik  (2L  -  6z)  , 


r  2  2  1 

B  =  (l/2)y[<5z  (a  +  6z  )  -  lj  . 

For  E*  perpendicular  to  the  cylinder,  the  diagonal  elements  are 


C  =  1  -p-h 
nn  nn 


(3-31) 


(3-32) 


where 


h  =  i  Ya  6z  e^/^a 
nn 

3.  2. 5  EXTERNAL  FIELDS 

T 

The  Reid  outside  a  cylinder  is  found  by  evaluating  Equation  3-1  with  the  values  of  E 
determined  by  solving  Equation  3-30.  The  integrations  are  the  volumes  of  the  discs  to 
preserve  dimensional  units  of  die  field,  not  on  set  of  points  corresponding  to  disc  cen¬ 
ters.  For  the  divergence  we  utilized  finite  differences. 

First  let  us  consider  the  contribution  of  the  vector  potential  to  the  external,  scattered 
field.  From  Equation  3-1, 


E®  (r»)  =  I.j/4* 


(3-33) 


where 


il  -  k 


j \k  -  1)  g  ET  dv  , 


(3-34) 


T  -1 

and  E  is  the  total  internal  field,  g  is  r  expikr  with  r  the  distance  between  integra¬ 
tion  and  observation  points,  and  k  is  the  dielectric  constant.  Equation  3-33  omits  the 
term 


E  =  V»X,/4jt 

“S 


(3-35) 


•  vwi 


(3-36) 


with  the  prime  indicating  differentiation  on  the  observation  coordinates  where 


k 


/ 


g  V  .  [(k  -  1)  E]  dv 


The  result  from  1^  is  for  the  q  component  (q  =  y  or  z) 


E®  <*!^  *’)  =-j(ka)2  (k-1)  e11*  E^  (z)  e  n  ^  [x'2  +  (z^-z')2] 


-1/2 

(3-37) 


For  I.  the  integrations  are  different  for  E*  and  E1 . 
2  y  z 


To  evaluate  I  for  E  ,  in  Equation  3-35  we  utilize  polar  coordinates.  So 
“2  y 
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=/‘[£<‘-1,epV;£(‘-1>e»\K 


(3-38) 


Because  the  dielectric  cylinder  has  small  diameter,  we  assume  3E p/dp  is  negligible. 
This  assumption  was  also  made  in  computing  the  internal  field.  Thus 


k  *  "(K 


-x,/ 


g  E^  6  (a)  p  dp  d$  dz  -  (k 


cos  <p  +  E  sin  <{>)  dp  d$  dz 
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(3-39) 


The  Dirac  delta  function  comes  from  differentiating  (k  -  1).  Now  the  incident  field  has 
only  the  El  component,  so  we  assume  E  is  negligible.  After  some  manipulation,  we 


have 


k  =  l2k{K 


-•>■77 

J-L  Jo 


2  it  k7h  1 

E  e  b"  cos  (tf>  -  <p')  sin  <p  <H>  dz  ,  (3-40) 
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2  2 

where  2L  is  cylinder  length  and  b  is  p'  +  (z  -  z1)  .  By  stationary  phase  we  integrate 


on  z  to  obtain 


E®  =  i|«c-  1)n/T (ka)2  cos  eikp'  (kp')3/2  (F+  +  F_)  Ey  (3-41) 


where  F  is  the  Fresnel  integral  evaluated  at  (2/Xp1)  (L  ±  z'). 


,.1/2 


For  Ez,  that  is  E  parallel  to  the  rod,  the  scattered  field  is  associated  with  delta 
functions,  at  the  ends  of  the  cylinder,  arising  from  the  divergence  operation.  The  re¬ 
sult  is  assuming  the  field  constant  in  the  end  cells 

s  1  2  ikD  -1  -2  -i 

E2(x»,  z1)  =  -EN(K-l)a  e  D  {LD  -  ikLD  1  ,  (3-42) 


where 


D  is  <p'2  +  L2) 

for  E1 , 

y 

ES  (x»,  z*)  =  i  (k  -  1)  l  s^r  (ka)2  cos  <?>*  eikp  (kpf  3/2  (F  +  F  )  E  ,  (3-43) 

y  o  +  ~  y 

if  we  assume  the  field  is  constant  in  tie  cylinder,  where  F±  is  the  Fresnel  integral 
evaluated  at  j2/\p ’  (L  ±  z'). 

3.3  COMPUTATIONAL  EXAMPLES 

The  formulas  of  Sections  3. 2. 4  and  3. 2. 5  were  utilized  to  calculate  the  fields  near  a 
cylinder.  The  first  step  was  to  invert  the  matrix  in  Equation  3-30  to  determine  the 
total  field  within  the  cylinder.  The  second  step  was  to  compute  the  external  field. 

Figure  3-2  shows  |  E  |  inside  a  dielectric  cylinder  for  N  =  26,  with  k  equal  2. 6,  di¬ 
ameter  0. 139",  length  1.  9".  The  wavelength  was  1. 259",  and  the  electric  field  was 
parallel  to  the  cylinder  axis.  The  figure  gives  |E^|  computed  according  to  the  theory 
of  Section  3. 2. 4  with  A  and  B  in  Equations  3-31  and  3-32  set  equal  to  zero  and  again 
for  A  and  B  evaluated  according  to  these  equations.  We  see  that  the  effect  of  A  and  B 
is  to  change  |Et|  near  the  ends  of  the  cylinder.  Figure  3-2  also  shows  the  exact  val¬ 
ue  of  |Et|  for  an  infinitely  long  cylinder  possessing  diameter  and  dielectric  constant 
equal  to  those  of  the  finite  cylinder. 

The  external  field  was  computed  first  from  the  vector  potential  formula.  Equation 
3-37  with  E^  five  internal  field  computed  from  both  the  vector  and  scalar  potentials. 
The  total,  external  field  was  computed  for  E  parallel  to  the  cylinder  axis  for  distance 
p'  equal  X/2.  See  Figure  3-3.  The  total,  external  field  for  E  perpendicular  to  the 
cylinder  axis  is  shown  in  Figure  3-4  again  for  p*  =  X/2.  The  effect  of  adding  the  con¬ 
tribution  from  the  gradient  of  the  scalar  potential  also  is  shown  in  Figures  3-3  and 
3-4. 
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Figure  3-2.  Total  Internal  Field  IeTI  inside  cylinder,  k  =  2.6,  length  1.9", 
diameter  0. 139"  calculated  by  moment  method.  Exact  solution  for 
infinite  cylinder  ( - )  \  =  1.259",  N  =  26. 

3.4  MEASUREMENT 

The  field  near  the  cylinder  of  Section  3-3  was  measured  with  a  dipole  probe  in  a  mi¬ 
crowave  interferometer.  Reference  1  described  apparatus  and  procedures.  Figures 
3-3  and  3-4  show  measured  data.  The  measurements  were  made  with  the  cylinder 
suspended  by  Nylon  threads.  Additional  data  were  given  in  Reference  1. 

3. 5  DISCUSSION 

In  general,  the  discrepancies  better-  «he  measured  and  computed  nearfield  values  of 
ET  are  reasonably  small  and  verify  the  theory,  the  mathematical  approximations,  and 
the  numerical  procedures.  Figure  3-3  shows  die  scalar  potential  contribution  negligi¬ 
bly  changes  the  oomputed  values  for  the  field  parallel  to  the  rod.  Discrepancies  are 
small  even  without  the  scalar  potential  contribution.  Figure  3-4  shows  that  the  scalar 
potential  contribution  is  significant  for  fits  electric  field  orthogonal  to  the  rod. 

The  experimental  data  show  |  E  ~  | 2  is  larger  for  file  field  parallel  to  the  cylinder  than 
for  file  field  orthogonal.  This  result  is  intuitive.  The  phase  values  near  the  rod  cen¬ 
ter  are  approximately  equal  for  both  oases.  Near  the  ends  of  the  cylinder  phase  de¬ 
creases  more  rapidly  for  the  parallel  oase. 
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Figure  3-3.  External  Total  field  for  1. 5  X  long  rod;  k  =  2. 6;  diameter  =  0. 1J 
X  *  1. 259  at  x  =  0.  5  X.  Measured  (x);  computed  from  vector  potential  (• ); 
from  vector  and  scalar  potentials  (+).  Electric  field  parallel  to  rod. 

The  dashed  lines  show  the  ends  of  the  cylinder. 


4.  NEARFIELD  SCATTERING  BY  A  HOLLOW  DIELECTRIC  CYLINDER 

This  section  describes  nearfield  scattering  by  hollow  dielectric  cylinder. 

4.1  THEORY 

The  integral  equation  for  the  electric  field  is  obtained  from  the  definition  that  the 
scattered  field  is 

5s  ■  -sri ■v'ti  <4-« 

where  the  vector  potential  is  a  volume  integral 

A  =  -  i wpQ  (4tt)  7  (k-1)  g  E  dV  ,  (4-2) 

the  scalar  potential  is 

*  -  -  <4irf  7  g  ?  •  (x-1)  E  dV  ,  (4-3) 

and  g  is  r_1  expikr,  with  r  the  distance  between  observation  and  integration  points  k  is 
the  dielectric  constant,  and  the  prime  denotes  differentiation  on  the  observation  point 
coordinates.  In  these  formulas  Eis  the  total  field,  which  is  unknown.  The  potentials 
arise  from  a  polarization  current 

Jpol  *  -  l"€o  I"'1*  5  (4-4> 

and  a  polarization  charge 

Ppol  =  -  V  .  (K-1)  E  .  (4-5) 

where  the  divergence  is  over  the  integration  point  coordinates. 

The  integral  equation  is 

E  -  ES  =  E1  (4-6) 

g 

where  E  is  given  by  Equation  (4-1)  with  A  and  $  from  Equations  (4-2)  and  (4-3). 


4-1 


The  integral  equation  is  solved  for  the  interior  of  the  dielectric  cylinder.  The  equa¬ 
tion  is  changed  to  a  set  of  simultaneous  algebraic  equations.  The  cylinder  is  divided 
into  rings  and  subdivided  into  cells  as  in  Figures  4-1,  4-2  and  4-3.  The  integral 
equation  is  evaluated  at  the  center  of  each  cell  to  generate  as  many  equations  as  there 
are  cells.  At  the  center  of  a  particular  cell,  labelled  with  index  m,  the  integral  equa¬ 
tion  for  each  rectangular  component  is 


g 

for  l<msN,  where  N  is  the  number  of  cells.  E  is  a  sum  of  contributions  from  all 

«  m 
cells,  so 


E  ~y  ES  =  E1 
m  La  mn  m 


Y 


iOHOOS 


Figure  4-1.  Dielectric  Ring  and  Coordinate  System 


This  equation  can  be  put  into  obvious  matrix  form  by  writing 


Figure  4-2.  Subdivision  of  Hollow  Cylinder  Into  Rings.  An  observation  point  has 
coordinates  (x' ,  y* ,  z' );  it  may  be  inside  or  outside  the  cylinder. 
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Figure  4-3.  Coordinates  for  a  Cell.  A  typical  cell  has  dimensions  as  follows: 
radial  26p,  circumferential  26<p,  and  longitudinal  26 z 
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where 


C  =  1  -  (E  /E  )  , 
mm  '  mm  m  / 

and 

C  =  -  E®  /E 
mn  mn  n 

g 

The  next  step  is  to  evaluate  the  scattered  fields  Emn.  Formulas  are  derived  in  the 
following  sections.  Diagonal  and  off-diagonal  matrix  elements  are  separately  given, 
and  the  vector  and  scalar  potentials  are  separately  evaluated. 

4.2  SCATTERED  FIELD  FOR  DIAGONAL  ELEMENTS 

4. 2. 1  VECTOR  POTENTIAL 

The  scattered  field  from  the  vector  potential  is 

ES  =  k2  (k-1)  (4 it)"1/  g  E  dV  .  (4-7) 

Let  us  assume  horizontal  polarization  both  incident  and  within  the  ring.  This  assump¬ 
tion  may  be  valid  if  we  consider  thin  rings  and  axial  incidence.  For  the  cell  labelled 
m  we  have 


E*  *  fc2  <47r>_1  f  8  E  dV  (4-8) 

xm  J  x  m 

To  simplify  notation  let 

IAn,  =  /8Ex'IVm  <4-9) 

Since  the  integration  extends  only  over  the  cell  labelled  m, 

f  2  2  21^ 

r  =  [p  +  Pm  -  2ppm  cos  (<p  -  <*>m)  +  (z  -  zm)  J  (4-10) 

with  integration  limits  pm  -  5p  £  p  s  +  6p,  zm  -  <5z  s  z  <  zm  <5z,  <pm  -  <50  i  0 
s  <pm  +  6<p ,  Introduce  abbreviations  0  =  0-  0m,  t  -  z  -  zm,  and  <r  *  p  -  pg v  Then 
on  assuming  0  -  <pm  is  small,  we  have 
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so  that 


bz  b<f>  bp  ikr 

f  f  f  E  r  e  m  (p  +  <r)  cos  0  dtr  d0  dt 
J  J  J  x  m  tn 

-5z  -6<p  -bp 


(4-12) 


The  factor  cos  0  in  the  integrand  arises  from  first  expressing  the  x-component  of  field 
in  polar  coordinates,  expanding  cos  (0  +  <f>m)  and  sin  (0  +  <pm),  and  then  using  small 
angle  approximations  for  sin  0  and  cos  0.  Now  assume  rm  is  small;  that  is,  e  and  £ 
are  small.  Moreover,  assume  Ex  is  constant  Em  in  the  cell.  Thus  we  can  expand  the 
exponential,  to  obtain 


I  =  I(1)  ♦  !<*>  ; 
Am  Am  Am  ’ 


(4-13) 


where 


6z  6<f>  bp 

^m  =  ikEm  J  /  /  (pm  +  *>  008  9  * d9  «  <4"14> 

-<5z  -6<p  -bp 


This  integral  is,  if  we  omit  a  in  p  +  <t, 

m 


1^  =  ik  (8  bz  bp  p  6<f>)  E 
Am  m  m 


The  factor  in  parentheses  is  the  volume  of  the  cell. 


(4-15) 


The  term  I;  '  is  more  complicated.  We  omit  limits,  which  equal  those  in  Equation 
(4-12)  Am 


■<l)  - 
Am 


Iff  fan  *  "»  "m  **  *  *  *  (p^  +  o)  cos  0  do-  d0  d£  E^  (4-1 


The  integral  over  £  leads  to  a  logarithmic  function  so 

I?*  *  f  f  (P  +  <r)  cos  9  log  [(R+l)/(R-l)j  do-  d0  E  (4-17) 

Am  J  J  m  m 

where  R  is  £l  +  (s*/^)]1^2.  and 
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where  s+  is  [<r2  +  (pm  +  o)  pm  d‘ 
root,  we  obtain 


fll/2  2  2 

j  .  Assuming  s+  <  5z  ,  to  expand  the  square 


JAm  =  1 1  (Pm  +  008  lo®  I45*2  s+2  [*  +  (s+/26z)2]}  do  d0  Em  (4-18) 

c  c  2 

-  (p  +  <t)  cos  0  log  (25/s  )  do  d0  E 
J  J  m  +  m 

+  //  (s+/26z)2  +  a)  cos  6  do  dd  E^ 

We  assume  that  cos  6  **  1  and  o  <k  pm  to  carry  out  the  integrals.  We  finally  have 
that  the  x  component  of  the  scattered  field  is 


I  =  it  a  [1  -  log  (s/25z)  ]  +  p  — ^  s  , 
Am  2  m  2  2 

36  z 


(4-19) 


where 


2  x  2  2 

s,>  =  fip  +  Pm6<*> 


Equations  (2-8),  (2-13),  (2-15),  and  (2-19)  give  the  scattered  field  in  the  cell  labelled 
m 


EAm  '  k2  "<-1) 


i  -  k  5z  5p  p  6<p  +  pm-- 
*  m  12x6 z2 


6p5(t>  s2  +  -  j  [l  -  log  (s2/26z)2j  J 


m 


(4-20) 


For  brevity  we  omit  the  details  of  the  integrations  over  o  and  0.  The  integral  was 
done  by  integrating  over  local  polar  coordinates.  In  effect  then  our  cells  are  small 
right  circular  cylinders  not  the  truncated  wedges  defined  by  pmd0do.  This  approxi¬ 
mation  procedure  is  described  near  the  end  of  Paragraph  4. 2. 2. 


4.2.2  SCALAR  POTENTIAL 


The  scattered  field  from  the  scalar  potential  is 


~s 


(4-21) 
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where 


(4-22) 


I2  -  /gV  •  (K-l)EdV  . 


To  represent  the  dielectric  constant  we  use  the  unit  step  functions  as  functions  of  p 
only.  That  is,  we  shall  understand, 

(t- 1)  =  (a- 1)  (H  (Pm  -  6p)  -  H  <Pm  -  5p)J  (4-23) 

where  on  the  right  side  (k-1)  is  a  constant,  but  on  the  left  k  is  a  function  of  p.  Be¬ 
cause  we  have  assumed  Ez  =  0,  it  is  unnecessary  to  specify  the  z  dependence.  Final¬ 
ly,  k  is  assumed  independent  of  <p.  Therefore 

ikr 

l2  =  J^T  V  *  (K"1)  ^  (Pm  '  6p>  "  H  <pm  +  6p>l  I  (I,,  cos  <f>  -  1^  sin  <f>)  dV 

(4-24) 


where  1  and  1  are  unit  vectors  in  the  r  and  <t>  directions  respectively. 

“T  “TO 

We  assume  div  E  =  0,  in  the  cell  so 

where 


(4-25) 


llT)  =  /  g<5  (p  -  <5p)  cos  <|>  dV  , 

J  m  m 

(4-26) 

42)  =  / (Pm  +  6P)  cos  0  dVjjj  . 

(4-27) 

=  J  g  cos  <p  dp  &t>  dz  . 

(4-28) 

r(l)  t(2) 


t(D 


The  first  two  integrals,  and  ,  approximately  cauicel.  To  see  this  consider  Ij>  . 
Assume  kr  is  small  and  expand  explkr,  Thus 


I21)  =  /  1  +  ik)  6  (p^  -  6 p)  cos  <f>  dV^ 


4-7 


■> 

atmmmm 


The  second  term  of  die  integraiid,  with  factor  ik,  gives  -4ik  (Pj^Sp)  <5z  6<f>  cos  $m. 
This  term  is  omitted  because  it  approximately  cancels  the  term  -4ik  (pm  +  6p)  6z  6 <p 
that  arises  from  v£\  The  first  term,  with  factor  r-1,  gives 

41}  “  (X  "^)  cos  /  [(£/pm)2  +  (fip/pm)2  +  (J  “^r)  cos  d«  dC  ; 

(4— 2£ 

recall  0  and  £  were  defined  earlier.  Following  integration  on  £  we  have 

n 

^1  - 'j  cos  J(fi nP+-inP  )  cos  0  d0  d£ 

(4-3C 

where 

.  (p  5p)  e ^ 

p  _  q  =fc  6z  (  vpm  p' 

±  2q 

xn 

(4-31 

and 

/  2  2 
q  =  a  6z  +  6p 

(4-32 

(21 

The  evaluation  of  I' 

proceeds  in  the  same  way.  The  result  is  that 

42)  - 

^  cos  y<enQ+  -in  Q  )  cos  0  d0  d£ 

(4-32 

where 

.  .  .  (p  +  6p)  0^ 

o  «  (q  *  Sz)  'pm 
*  Pm 

(4-34 

Now  consider 

ll1)  -  l£2)  -  f  in  (L./L  )  d0 

(4-39 

where 


s2  -  6z2  „2  'm 

- —  +  0  + 


2  A  2 


fl4±3zj5p  q2 


m 


4s 


8  P, 


(4-36) 
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Inspection  of  L+/L_  shows  that  it  tends  to  unit  value  for  small  0.  For  example,  at 
0=0  the  value  is  1;  for  0  =  0. 2,  the  value  is  1. 046.  Therefore  the  logarithm  is 
small  so  the  contributions  of  and  I^2)  from  are  omitted. 

(3) 

To  evaluate  ,  expand  the  exponential,  expikr,  in  the  integrand  as  1  +  ikr.  The 
first  term  is  omitted  because  it  is  independent  of  the  observation  point  coordinates 
so  the  V  operator  yields  zero.  From  the  second  term  we  obtain 


r<3) 


h'1 


cos  <p  dp  <ty>  dz 


(4-37) 


Again  we  use  substitution  z-z  =  £,  p  -  p  =  a,  <f>  -  <f>  =  0  to  obtain 

m  m  m 

-1/2 

Ig3^  =  cos  j  J  j  [ff2  +  (Pm  +  a>  Pm  6 2  +  ^2]  008  0  dp  d0  d£  .  (4-38) 

We  assume  0  is  small  so  cos  0=1.  Integration  on  £  gives 


43>  "  008  *m 


jyj 

J  J  lL  4  6  z  J  (p^  +  0)  pm  0  + 


m  m 


da  d0 


(4-39) 


Because  wall  thickness  is  much  smaller  than  radius,  we  omit  the  term 


a  p  0  /46  z 
m 


in  file  integrands.  The  integrations  on  a  and  0  give 


43>  *  ' “m'1  “*  *m  [’  '  l0*  <"m/24  z>2] 


-2 

2  6p  6$  (6*)  2 

s  + - —z -  s  oos  p 

m  3  m  m 


(4-40) 
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To  obtain  the  field  apply  the  V'  operator  to  1^  .  The  p  and  <p  components  are  as 
follows: 

(EL)p  -  -  ^l>  j5!2  [*  +  (*  -  (!)2  5*'2)  ( 1  - T  "m'2  S*'2) 
6p  pm  **  1  ) 

+ - T  cos  <p  E  i 

3,  J  m  *m 


(EU  ■-**>(<)*  ( 


2-i  p  <5p  6<p 

_  +  JE - — -  sin  0  E  ( 

*■*  3rf*2  I  m  “ 


From  Equations  (4-41)  and  (4-42),  the  contribution  of  the  scalar  potential  to  the  x 
component  of  the  scattered  field  is 


EL  (*)  =  -  (K-l)  ^ 

xxn  2 


f1*^1 1- i)*f] 


-4—)  +  fj  COS2 

2p  6<p  >  J 

m 


where 


*{wj  <d  +  *>  sln'  V  Pan 


d  =  1 


- Io*  (f ) 


f  ■  p  fip  66/(3 k  6z) 


4.3  SCATTERED  FIELD  FOR  OFF-DIAGONAL  ELEMENTS 
4. 3. 1  VECTOR  POTENTIAL 

The  scattered  field  at  the  cell  centers  results  from  a  finite  sum  of  integrals  over  cells 
like  those  in  Figure  2-3.  The  vector  potential  gives  for  the  x  component 

E®  -  k2  (4jt)“1  (k-I)S’  E  I  ,  (4-46) 

Am  n 


where 


I  =  f  g  (1  cos  0-1  sin  <p )  d0 

■I  p 


(4-47) 


and  the  prime  suggests  the  sum  excludes  the  term  for  n  =  m. 


In  g,  set  p  -  p  ;  for  small  0, 
n  m 


=  f  2b  p2  +  (z  - 
[  "m  n 


z  )  +2p  bo  +  2  (z  -  z  )  £ 
mm  n  m 


j!/2 


(4-48) 


where  b  is  1  -  cos  (<pn  -  4>m)-  The  binomial  expansion  gives  approximately  that 

-1 


r  =  /3  +  p  (p  bo-  +  (z  -  z  )  £]  , 
m  n  m 


(4-49) 


where 


f  2  2l1/2 

=  [2bPm  +  <zn  *  zm>  J 


With  Equation  (2-49), 


-i  r  lk  <zn  "  zra>  c  /P 

1*0  e  fe  e  (p  +  o)dod0dt. 

j  n 


(4-50) 


(4-51) 


8  p  6P6Z600"1  elkp 
m 


/  \2  2 
1  -  (1/8)  (k  Cm  M  flS'1)  -J  [k  <*„  -  V  41  *  J 
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The  higher  terms  have  magnitudes  much  less  ft  an  l,  so 


I 


-1  lkfi 

8  6p  6z  5 ip  P  e 


(4-52) 


The  scattered  field  is 


s  2  2  -l  ikfl 

EAm  =  k  P«  5*  5z  6*  2 ' ’  Exn  P  6 


(4-53) 


4.  3. 2  SCALAK  POTENTIAL 


The  derivation  of  matrix  elements  is  slightly  more  general  for  this  case  because  of 

the  gradient  on  observation  point  coordinates.  That  is,  the  gradient  operation  V*  is 

performed,  after  that  set  p  =  p  .  The  distance  in  the  expression  for  g  is 

m  n 


-  [<■  *  fj  -  2» 


p  cos  (<p  -  <t>  )  + 
m  m 


i\1/2 

<‘-VJ  • 


Now  let  z  =  2  +  f ,  p  -  p  +  a  so  that 
n  n 


r  = 


[s2  +  2p^  a  -  2p 


1 


1/2 


ff  cos  (<#>  -  4>  )  +  2  (Z  -  Z  ) 
m  n  m  n  mj 


(4-54) 


(4-55) 


where 


'  (f 


2  2  2  ' 

S  ■  lpn  +  Pm  +Znm  ~2pnPm<T 
n  m  run  n  m  j 


1/2 


The  binomial  expansion  gives  approximately  that 


(4-56) 


r  -  s  +  p  -  p  cos  (<p  -  <f>  )  +  z  s  t  . 
n  m  n  m  nm 

Now  consider  of  Equation  (4-28).  With  r  from  Equation  (4-57), 

1^  ■  8  6p  dtp  6z  eikS  s  1  cos  <p 


(4-57) 


(4-58) 


if  we  approximate  the  sine  functions  that  result  from  integration  by  unit  value. 


From  Equations  (4-26)  and  (4-27), 


T(U  t(2) 

2  "  2 


8  6z  &<(>  6p  e*  s  * 

(3) 


1  +  ikp  [p  -  p  cos  <f>  -  <#>  )]  s 
n  n  m  n  m 


-1 


.(4-59) 


The  scattered  field  from  1^  for  the  nth  cell  is 

E8^  =  -i^E  V’I<3> 
n  2  4jt  xn  2 


=  (k-1)  6p  6<t>  6z  cos  <p  E  V*  eikS  s  1  E 

7  n  n  xn 


(4-60) 


The  differentiations  are  on  p^  and  <p^.  The  x  component  is 


E"  1^  =  -—  (x-1)  <5p  6z  6<f>  cos  <p  (ik  -  s  1)  eiks  p  s  2  sin  «>  -  <p  )  sin  <f>  E 
x  &  7r  n  m  n  m  m  xn 


The  gradient  operation  on  1^  -  1^,  gives  two  terms;  in  fact. 


(4-61) 


(4-62) 


where 


A  E  = 


iks 

s 


[■ 


2  -2  ,  2 
s  s  2p  s  b-l-ip 
nm  m  m  m 


bks  *  I  +  -5?  cos  <p  b  T2  p  bs  * 
J  s  m  L  m 


+  (1  -  b)  -  ik  pmb] 


(4-63) 


with 


s 

nm 


sin  (^  -  <p  )  and  s 
n  m  m 


sin  4 

m 


4.4  EXTERNAL  FIELDS 

The  field  outside  the  dielectric  region  is  found  by  evaluating  Equations  (4-7)  and 
(4-21).  The  integration  extends  over  the  dielectric  region.  It  can  be  done  after  the 
total  field  inside  the  dielectric  is  determined  by  solving  the  integral  equation. 
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The  field  near  a  cylinder  is  evaluated  to  test  the  approximations  in  the  derivations. 

The  near  field  also  is  useful  in  understanding  the  effects  of  a  dielectric  cylinder  or  a 
radome  on  an  antenna. 

The  field  from  the  vector  potential  is 

ES  (A)  =  i  cd  A  . 

A  measurement  requires  a  probe.  We  consider  a  half-wave  dipole  and  weight  the  con¬ 
tribution  of  each  cell  by  the  far-field  pattern.  To  simplify  the  expressions,  consider 
the  field  on  the  z-axis.  The  sector  potential  gives  for  the  x  component  (from  one  ring) 

ES  (A^)  =  (k-1)  k2  (4n)  1  Jg  E^  cos  cos  0)  (sin  0)  1  pdpd  <t>dz  (4-64) 

where  the  factor  cos  (r/2  cos  0)/sin  0  is  the  far-field  pattern  of  the  probe.  Consider 
the  distance  large  enough  so  that 


r  =  R  +  (an  -  a p  cos  <p  -  zzT)  R  1  (4-65) 

1/2 

where  a  is  the  mean  radius  of  the  cylinder  and  R  is  (a2  +  z'2)  .  Furthermore  ap¬ 

proximate  sin  0; 


2  2 

sin  0  =  1-  (a/R)  cos  <£>  (4-66) 

To  develop  a  simple  formula  we  omit  the  dependence  of  E  on  <p  and  use  an  average 
value  Eq.  The  integrations  are  elementary,  so  that 

s  2  ikR  -1 

E  (AJ)  ■  2  (k-1)  k  a  6p  <5z  Eq  e  R  sine  (ka  6p/R)  sine  (kzf  6z/R)  P 


[l  -  (a/2R)2  (j-  -  l)j  . 

P  is  [1  -  (a/2R)2  (tt2/4  -  1)],  the  probe  correction. 

The  field  from  the  scalar  potential  is 

E®  (*)  =  (47T)"1  7f  fgV  .  (K-1)  E  dV  . 

X  *  X 


(4-67) 


(4-68) 
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The  procedures  are  like  those  in  Section  4. 2. 2  and  4. 3. 2.  So 


£»<*) 


2  ikR  , 

-K^DfxH  E0  /e“<*")  f--.]  cos%  dcrd*  dz  ,  (4-69) 

R 


where 


(...)  =  (ao  -  a p'  cos  <p  -  zz')  R 


(4-70) 


and 


-1. 


I. . .  ]  =  [5  (cr  +  6p)  -  5  (a  -  6p)  -  a  ]  . 
After  doing  the  integrals,  we  have 

-1 


(4-71) 


K  <*i>  ’ 1 ’  w 6p  “im  (r2x)  (I + x)  Jo  (x) (nr)  p'  E< 


(4-72) 


2  2  2 

Note  the  probe  correction  Pf  is  [1  -  3/8  (ir  /4  -  1)  a  /R  J;  it  difFers  from  that  in 
Equation  (4-67). 

The  total  field  is 


E^  =  E1  +  E®  (A)  +  E®  (*) 

a  A  A 

T  2  T 

The  observable  quantities  are  intensity  |  E^  |  and  phase,  the  argument  of  E^ . 


(4-73) 


This  theory  does  not  include  the  effects  of  fields  reflected  from  the  probe  to  the  scat- 
terer  and  scattered  back  to  the  probe. 


The  formulas  generalize  readily  to  several  rings  forming  a  cylinder.  For  example 
for  two  rings 


s  .  „s  ..  .  /_  ikz'  6z/R  _  _ikz*  6z  1R 


e;<v  * E"  (v  (E»i  • 


+  e 
02 


) 


(4-74) 


where  Eg^  and  Eq2  are  the  fields  in  each  of  the  rings  and  E  (Aj)  is  from  Equation 
(4-67).  The  scalar  potential  gives 


The  first  example  was  for  a  single  ring  diameter  2a  equal  X,  length  2  6z  =  0. 16  X,  di¬ 
electric  constant  2.6.  Because  incidence  was  axial,  necessary  conditions  on  the  com¬ 
puted  results  were  that  the  field  at  <p  ~  0  equal  that  a.t<p  -  r  and  that  the  field  at  <p 
»  7r/2  equal  that  at  <f>  =  3  ir/2.  This  condition  was  well  satisfied.  Figure  4-4  shows 
computed  internal  values  of  ET  for  several  values  of  N. 


Figure  4-5  shows  computed  internal  values  of  ET  for  a  cylinder,  diameter  X,  length 
0. 32  X.  Two  rings  were  utilized. 

T 

Figure  4-6  shows  computed  values  of  E  outside  the  ring,  on  the  z-axis,  with  and 
without  probe  corrections.  Figure  4-7  shows  computed  values  of  ET  outside  the  cyl¬ 
inder.  Figure  4-8  shows  computed  values  for  the  cylinder,  but  the  area  of  the  cylin¬ 
der  cross  section  was  reduced  by  80%  in  evaluating  the  external  field  integrals. 


T 

Figure  4-4.  Phase  and  Amplitude  of  E  for  Single  Ring:  6z  *  0. 100  in. ,  6p  ■  0. 030 
in. ,  x  *  2.6,  frequency  ■  9.375  GHz.  For  $  ■  0*  (o),  <p  ■  90*  (•) 
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Figure  4-6.  Total  Field  Behind  Dielectric  Ring.  Length:  0. 2  in. ;  thickness  26  p: 

0.06  in.;  mean  diameter  2a:  1. 25  in.;  dielectric  constant  2.6.  Meas¬ 
ured  with  dipole  probe:  (•);  computed  without  probe  corrections  (o); 
computed  with  probe  correction  (+).  Internal  field:  1.19  exp  i  7. 8* . 
tor  N  *  64 
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